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Abstract: A new approach to Gamma/Hadron separation algorithms is proposed. The differences be- 
tween Gamma and Hadron showers are notorious in two main aspects. The first is the wideness of the 
shower, and the second is the distribution of the angles of emission of Cherenkov photons in respect to 
the shower main axis. Using more than one IAC telescope, and their respective bi-dimensional images 
of arrival directions of the Cherenkov photons, the 3D geometrical characteristics of the shower can be 
reconstructed. 



Introduction 

In recent years, the search for Gamma Ray sources 
has become one of the main concerns for astro- 
physics. The development of ground-based Imag- 
ing Atmospheric Cherenkov Telescopes (IACTs) 
allowed the detection of Gamma Ray showers in 
the Very High Energy domain (in the order of the 
TeV) [1, 3]. This opened a possibility for new 
physics to be studied. 

The number of known emitters has greatly in- 
creased in the last years. However, most of them 
remain unidentified and the prospect of discover- 
ing new objects is good. Several telescopes were 
developed to study these sources and, simultane- 
ously, search for new ones. 

Most of the IACTs already built, as well as the ones 
planned for the near future, are characterized by 
having a Field of View (FoV) of only a few de- 
grees. While providing good resolutions for the 
study of the known sources, a sky survey of a con- 
siderable portion of the sky would consume too 
much time. 

A new approach to the IAC technique is being pro- 
posed by GAW [2], that will introduce a large FoV 
up to 24°x 24°. One of the main goals in GAW 
is the serendipity search of sources, by performing 
a sky survey of 360°x 60°. In order to effectively 
locate sources, a good shower reconstruction algo- 
rithm has to be used. Furthermore, it is also neces- 



sary to distinguish the 7-ray showers from the ones 
induced by other primaries. 

In GAW, three telescopes will be placed at the ver- 
tices of an aproximately equilateral triangle of 80 
m side. This configuration was used as a case study 
to test the proposed method. The origin of the ref- 
erence system is located at the centre of the trian- 
gle. 

Generally, a bi-dimensional analysis of the images 
obtained by the IACTs is performed to calculate 
the Hillas parameters [3]. These are then used to 
characterize the detected shower and determine the 
type of primary that originated it. The Hillas pa- 
rameters analysis is based solely on the projections 
of the shower in each telescope. Therefore, there 
is the possibility that additional information can be 
used by reconstructing the shower in 3D. 

With the positions of the telescopes and the ar- 
rival directions of the Cherenkov photons emitted 
in the cascade, a geometrical representation of the 
shower can be constructed. Different primaries 
have different characteristics, such as the shower 
being wider for proton induced showers than 7 in- 
duced ones. These differences can be translated 
in variables such as the impact parameter of each 
Cherenkov photon arrival direction (b), which is 
the minimum distance between the arrival direc- 
tion and the shower main axis. The shower wide- 
ness also relates to differences in the distribution of 
a, defined as the angle between the arriving pho- 
ton direction and the main axis of the shower. A 
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geometrical representation of both variables can be 
seen in figure 1 . 




Figure 1 : Reconstruction variables 

With these variables, it may be possible to both 
develop an iterative 3D method to geometrically 
reconstruct Extensive Air Showers seen by IACTs 
and effectively measure the difference between pri- 
maries to identify the primary of each event. 

The separation algorithm uses the distributions of b 
and a that are built with the reconstructed shower. 
The method proposed here is a preliminary one, 
that is still under test with simulated data. 

Geometrical shower reconstruction 

Geometrically, the primary particle path is a 
straight line defined by the direction and a position 
in space. The point chosen is located at ground 
level, being the point where the primary would 
hit if it crossed the atmosphere with no interac- 
tion, called the shower core position. The direction 
is characterized by the zenith (8) and the azimuth 
($). 

The developed reconstruction method of the geom- 
etry of the EAS uses the a, b variables discussed 
before. But in order to calculate them, it is neces- 
sary to have an approximate value of the core po- 
sition and the main axis. To solve this problem, 
an approximation is originally calculated in a first 
stage. On the second stage, the initial guess is sig- 
nificantly improved. 

As a charged particle moves through the atmo- 
sphere at a speed higher than the speed of light, 
it emits Cherenkov photons in a cone. It can be 



assumed that some of the photons that reached dif- 
ferent detectors were emitted by the same particle. 
With a set of N detectors (N> l),the intersection of 
photon directions from different detectors provides 
a guess for points belonging to the shower. In re- 
ality, the directions usually do not intersect and the 
closest point to both is used. 

The group of points obtained this way, as approxi- 
mated as the algorithm is, provide useful informa- 
tion about the shower as they closely follow the 
actual main direction of the shower. A first value 
can then be computed calculating the main iner- 
tia axis. Furthermore, they should be positioned 
around the axis, and their barycentre should belong 
to it. The initial core position is simply the point at 
ground level, that belongs to the line defined by the 
computed direction and the barycentre of the set of 
points. 

The second stage is an iterative process. In each 
iteration, the impact parameter is calculated for ev- 
ery photon. The sign of b is related to the position 
of the emission point, in respect to the plane de- 
fined by the shower axis and the detector the pho- 
ton belongs to. For each telescope, the average and 
standard deviation of 6s are computed. If the re- 
constructed shower line is close to the real one, the 
averages should be close to zero. Here, a problem 
arises with the photons that reach the detector with 
directions contrary to the main axis which results 
in b being calculated below the ground. Since this 
is not a valid solution, these points will not be used. 

At the end of each iteration, the set of photons used 
in each detector is reduced to the ones that are close 
to the average, within a few sigma. This guarantees 
that in the next iteration, the new axis and core po- 
sition are calculated from a group of points with 
less fluctuations. The method is iterated until the 
core position is stable enough (moving less than a 
few centimetres). 

To test this method, a hundred 7 showers for the 
energies of 800 GeV, 1000 GeV, 1500 GeV, 2000 
GeV and 3000 GeV were generated. The core po- 
sition was selected randomly from a square of 120 
m side, while 9 ranged from 0° to 30° and $ from 
0° to 360°. To approach reality, an efficiency of 
10% was imposed at the telescopes. The average 
of the absolute value of the errors in the determina- 
tion of the core position and the shower main axis, 
are represented in figures 2 and 3, respectively. 
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Figure 2: Errors for the position 




Figure 3: Error for the deviation angle in the direc- 
tion 



70 m and 150 m that represent small, medium and 
large distances. 

The interpretation of the a distribution is not 
straightforward. It reflects not only the emission 
angle but, since the air refraction index is a func- 
tion of height, also the spatial distribution of the 
shower electrons as well as the distribution of their 
momentum. In any case, we expect to have small a 
angles both for 7 and protons (mostly up to 1.5°) 
but wider distributions for protons as the protons 
showers are wider. 



Cherenkov Probability for 3000 GeV 
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The comparison of these results with simulated 
proton showers indicates that the method provides 
more accurate results for 7 induced showers than 
proton induced ones. 

Gamma/Hadron separation 

To analyse the differences, in terms of the pro- 
posed variables, between 7 induced showers and 
proton induced ones, one event of each was gen- 
erated. The energy of the primary particle was 3 
TeV for both cases, with a vertical direction. All of 
the Cherenkov photons emitted were used in order 
to build statistics. This corresponds to having an 
enormous detector on the ground with 100% effi- 
ciency. 

The study of these variables was performed by 
analysing the distribution curves for each variable 
at fixed radius. For each distance R, the photons 
that arrived inside the corona of R ± 0.5 m were 
considered. The radius presented here are 10 m, 
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Cherenkov Probability for 3000 GeV 
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Figure 4: a distribution for R = 10m, R = 70m, 
and R = 150m 
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For a, the distribution curves are shown in figure 
4. It is easily noticed that the curves are very dis- 
tinct at 10 m, overlapping at around 70 m, and only 
slightly separated at 150 m. This suggests that a 
can be used to separate the two types of showers 
when the shower is close to the telescope. 



Cherenkov Probability for Proton of 3000 GeV 
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Cherenkov Probability for Proton of 3000 GeV 
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point below ground. Then, wider showers will also 
be characterized by larger peaks at b = R. Fig- 
ure 5 represents the b distributions for the selected 
distances. 

The analysis of the image shows that the chosen 
solution for the points calculated below ground in- 
troduces a "pile-up" at b = R. For small dis- 
tances, the curves are almost indistinguishable as 
the "pile-up" dominates the distributions. How- 
ever, for medium and high values of R the "pile- 
up" becomes very useful. In these cases, while 7- 
ray induced showers have almost no photons there, 
the proton spectrum has a well-defined peak in that 
region. Therefore, it seems that b is a good variable 
to distinguish the two types of showers for medium 
and high distances. 

Currently, the good estimator that combines all the 
considered telescopes and both variables is under 
construction. Clearly, a dependence in R has to 
be introduced in order to select the good variable 
according to the distance as each has its own region 
of influence. Although not yet complete, this path 
of study seems to provide a good discrimination of 
7-ray induced showers. 

Conclusion and prospects 

The reconstruction method uses two spatial vari- 
ables that can be understood as the emission angle 
and the point of emission of each Cherenkov pho- 
ton. The tests performed so far indicate errors of 
about 10 m for the core position and a deviation of 
about 0.08° in the shower direction. 

The proposed variables for the separation of 7 and 
proton showers seem to allow an efficient separa- 
tion for the distances considered. An estimator is 
being built to account for the number of telescopes 
in use, the distance of each telescope to the core 
position and the good region for each variable. 



Figure 5: Distribution of b for R = 10m, R = 70m 
and R = 150m 

In what concerns b, its interpretation is simpler. It 
should reflect mainly the spatial distribution of the 
shower electrons and so, wider distributions are ex- 
pected for protons. One should note that b is set to 
R whenever the minimum b value is calculated at a 
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